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Abstract 

We examine here the nature of the central singularity forming in the spherically 
symmetric collapse of a dust cloud and it is shown that this is always a strong cur- 
vature singularity where gravitational tidal forces diverge powerfully. An important 
consequence is that the nature of the naked singularity forming in the dust collapse 
turns out to be stable against the perturbations in the initial data from which the 
collapse commences. 
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The final outcome of the gravitational collapse of a massive matter cloud is an issue of 

great interest from the perspective of gravitation theory as well as its possible astrophysical 
implications. When there is a continual collapse without any final equilibrium for the cloud, 
a black hole may form where the superdcnsc regions of the matter are hidden away from the 
outside observer within an event horizon of gravity or, depending on the nature of the initial 
data and the possible evolutions, a naked singularity could result as the end product of such 
a collapse. 

In this context, spherical dust collapse in general relativity has been studied quite ex- 
tensively. The prescription of matter as pressureless dust could perhaps be regarded as 
somewhat idealized. However, some authors have considered dust as a good approximation 
of the form of matter in the final stages of collapse^. All the same, the study of dust models 
over the decades has led not only to important new advances and insights into various as- 
pects of gravitational collapse theory, but has also laid the foundation of black hole physics. 
The class of solutions representing an inhomogeneous spherically symmetric dust cloud was 
discovered and studied by Tolman, Bondi, and Lemaitre^. In addition to its wide use in 
cosmology to model the universes which admit inhomogeneities, these models have been 
used extensively to investigate gravitational collapse in general relativity. For example, the 
special case of homogeneous dust ball was studied by Oppenheimer and Snyder"^, which led 
to the concept and theory of black holes. This model also inspired the cosmic censorship 
conjecture of Penrose^, which basically states that singularities which are visible to outside 
observers, either locally or globally, cannot develop in gravitational collapse from regular 
initial data in a stable manner. 

In the case of formation of a naked singularity in gravitational collapse, two aspects which 
become most significant are the strength of such a singularity in terms of the behavior of 
the gravitational tidal forces in its vicinity and its stability properties. Numerous papers on 
the strength have led to the classification of various subclasses within the Tolman-Bondi- 
Lemaitre (TBL) spacetimes, where the central singularity is either shown to be weak or 
strong. Thus there is considerable debate on the genuineness and physicality of the central 
singularity in dust collapse, especially for the cases where existing literature points towards 
the singularity being gravitationally weak (see, e.g., Refs [5-9], and references therein). The 
importance of this lies in the fact that even if a naked singularity occurs rather than a 
black hole, if it is gravitationally weak in some suitable sense, it may not have any physical 
implications and it may perhaps be removable by extending the spacetime through the 
same^°. Such, for example, would be the case for the shell-crossing naked singularities 
caused by the crossings of dust shells, through which a continuous extension of the metric has 
been shown to exist in certain cases, and the field equations hold in a distributional sense. 
Similarly, if a naked singularity is not stable in some well-defined sense, it may not have any 
significant physical consequences. The purpose of this paper is to investigate these aspects 
related to the central naked singularity forming in the dust collapse. As opposed to the 
shellcrossings mentioned above, the physical area of the dust shells vanish in this case when 
the matter shells collapse into the singularity at the center. We show here that whenever 
such a singularity forms as the end product of the collapse, it is always gravitationally strong. 
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Certain stability aspects related to naked singularities are then pointed out. 

The gravitational strength of a spacetime singularity is conveniently characterized in 
terms of the behavior of the linearly independent Jacobi fields along the timelike or null 
geodesies which terminate at the singularity in the past or future^^. In particular, a causal 
geodesic 7(5), incomplete at the affine parameter value s = Sq, is said to terminate in a strong 
curvature singularity at s = Sq if the volume three-form V{s) = Zi{s) A ^2(5) A Zs{s) (in the 
case of a null geodesic, this will be defined as a two-form) vanishes in the hmit as s — > sq 
for all linearly independent vorticity free Jacobi fields ^1,^2,^3 along 7(5). A sufficient 
condition for the causal geodesic 7(5) to terminate in a strong curvature singularity^^ is that 
in the limit of approach to the singularity, we must have along 7, 

(so - sfRabV^V' = const > 0, (1) 

where is the tangent vector to the geodesic. This expresses the requirement on the 
strength in terms of the rate of divergence of the Ricci tensor along the given trajectory. 
Essentially, the idea captured here is that in the limit of approach to such a singularity, the 
physical objects get crushed to a zero size, and so the idea of extension of spacetime through 
it would not make sense, characterizing this to be a genuine spacetime singularity (we refer 
to Ref. [12] for a further discussion). 

We now discuss the structure and strength of the singularity forming in the TBL collapse, 
an issue which is all the more important when it is naked. These models are fully character- 
ized by the initial data, which corresponds to two free functions, namely, the initial density 
function p(r) [or equivalently the mass function F{r)] and the velocity or energy function 
/(r), specified at the initial spacelike surface t — ti at the onset of collapse. Consider the 
dust collapse described by the TBL metric, 

ds^ = -dt'' + ^dr^ + R'dn^ (2) 

The energy-momentum tensor is that of dust 

= 6^15^,6 = e(t,r) = -^, (3) 

where e is the energy density, and R = R{t, r) is given by 

R' = fir) + (4) 

Here the dot and the prime denote partial derivatives with respect to the parameters t and 
r, respectively, and for collapse we have R < 0. Using the remaining coordinate freedom left 
in the choice of rescaling of the coordinate r we set at the initial time t — 0, 

R{0,r)=r. (5) 
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The functions F{r) and /(r) are the mass and the energy functions respectively, referred to 
above. We have 



p{r) = e(0, r) 




F(r) — / p{r)r^dr. 



(6) 



Since at the onset of collapse at the initial surface t = the spacetime should be singularity 
free, it follows that F{r) = r^h{r) with F{0) = 0, and that the initial density p{r) must at 
least be a function including the center at r = 0, i.e. 



where g{r) is at least C° function with g(0) — 1, n > 1, and pi > are constants. The 
term r" is the first nonvanishing higher order term in density at r = 0. Thus, if p'(0) ^ 0, 
then n = 1. For physical reasons we take that p(0) = pc ^ 0. The shell-focusing singularity 
appears at i? = along the curve t = ts{r) such that R[ts{r),r] = 0. The central singularity 
at r = appears (for the case / = 0) at the time 



It is known^'*^ that given any initial density profile for such a cloud, one can always choose 
a corresponding velocity function (or vice versa) describing the infall of the matter shells, 
such that either a black hole or a central naked singularity would develop as desired as 
the outcome of the collapse. Thus a naked singularity would occur in a wide range of 
TBL spacetimes depending on the choice of the initial data. The structure of this naked 
singularity for the marginally bound, time symmetric case, when the density was taken to 
be an even smooth function [i.e., p'(0) = and all the other odd derivatives of density also 
vanish] was analyzed analytically by Christodoulu^^, and later generalized to a wider class of 
TBL spacetimes with the same assumption of even smooth functions by Newman^, who also 
showed the singularity to be gravitationally weak along the radial null geodesies terminating 
at the central singularity. The formation and structure of this naked singularity was then 
analyzed^ with the only assumption of differentiability of the initial mass and velocity 
functions. Several later papers [see, e.g., Ref. [9], and references therein) analyzed the 
precise role of the first three derivatives at the center of the initial density function towards 
determining the nature and structure of this naked singularity. The conclusion prevailing 
presently from all of this analysis is that in cases when the first two derivatives of the density 
at the center were nonvanishing, the naked singularity was gravitationally weak, because the 
curvature strength along the outgoing radial null geodesies did not diverge sufficiently fast. 
In other words, the central singularity formed in a dust collapse is naked but gravitationally 
weak in cases (i) p'(0) < 0, (ii) p'(0) = 0, p"(0) < 0, however, gravitationally strong and 
naked in the cases when (ni) p'(0) = p"{0) — 0, with p"'{0) < and less than a certain 
maximum. 



p(r) = pc[l- pir''g{r)], 



(7) 
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(8) 
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Physically it is argued sometimes that the initial density distribution must be an even 
smooth function of r, therefore for all such cases the singularity would be weak. Furthermore, 
from the point of view of stability, any slight perturbation of a strong curvature singularity 
would involve the introduction of the term where p'{0) 7^ 0, p"(0) 7^ 0, and the formation of a 
strong curvature naked singularity would not be stable in this sense. As emphasized earlier, 
one could possibly remove a weak naked singularity by possibly extending the spacetime, 
which, however, is not possible for a strong curvature singularity. Though it has been shown^ 
that for any given arbitrary initial density distribution one could always choose the velocity 
function to be such that the resulting singularity is a strong curvature naked singularity, it 
is possible that the measure of such profiles may be vanishingly small in the space of all 
initial data, and in this sense such a strong naked singularity may not be stable. Thus it 
could be argued that from the point of view of both the stability and curvature strength the 
naked singularities forming in the dust collapse may have very limited bearing on cosmic 
censorship hypothesis. 

We shall show, however, that the structure of the central shell-focusing singularity is 
rather more complex than thought earlier. Considering the marginally bound case (/ = 0) 
for the sake of clarity, even the case of radial null geodesies can be seen to be quite involved. 
The geodesic equation in this case, using Eq.(2) and an integration of Eq.(4), is written as 



It is seen that in the limit of approach to the singularity both the numerator and denominator 
vanish in the above equation, and therefore even in the simplest case of radial null geodesies 
the nature of the singularity turns out to be rather complex, which is a node for the above first 
order equation, giving rise to a complicated topology of integral curves near the singularity. 
All the analysis and related conclusions so far on this naked singularity have been based on 
trying to understand only some of the famihes of these null geodesies. Because of the complex 
and difficult nature of the geodesic equations, especially in the timelike case, the discussion 
on them has been avoided so far in the literature. It is, therefore, proposed to investigate 
here the nature of timelike geodesies terminating in the naked central singularity. We find 
that this changes the existing perception on the nature of the central naked singularity as 
discussed above, and in fact we get timelike radial geodesies along which the curvature growth 
is powerful enough near the singularity to satisfy the strong curvature condition above. 

Let us consider the timelike causal curves in the TBL spacetime. Let U"- = {dx°'/ds) {s 
being the proper time along the trajectory) be the tangent to a timelike geodesic, satisfying 
C/"C/(j = —1. We can express radial timehke geodesies as 




(9) 



dr ^/F[r3/2 - |^t]i/3- 




(10) 
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where ± represent outgoing and ingoing curves and we have used the variable R instead of t 
in the last equation. The solutions to the set of two differential equations above in the form 
of R{t,r) ~ fi{r) and C/* = U*{R,r), describe the trajectories of particles following timelike 
geodesies. 

Because of the nonlinearity, solutions to the above differential equations are not available 
in general. However, there is an exact solution of the above described by 

= — = r=0. (12) 

This is an ingoing radial timelike geodesic 7(5) which in (t, r) plane is described by r = 0, 
terminating at the central singularity at the coordinate time t = ts{0), corresponding to a 
particle at the center r — following a timelike geodesic and terminating in the singularity 
in future. The tangent to this timelike geodesic is given by the above equation, where s is 
the proper time. The equation of the trajectory is simple and is given by 

ts — t — So — s, r — (13) 

where sq is the proper time when the particle crashes into the central singularity at r = 
at the coordinate time ^^(0). 

It is now possible to show that the above is not the only geodesic terminating in the 
future into the singularity. In fact, there are infinitely many families of future directed null 
and timelike geodesies from the past of the central singularity at r — 0,t — tg, terminating 
at the singularity in future. To see this, consider the past I^i'j) of the timelike geodesic 
7(5) given above. Since 7(5) is future endless and future incomplete, terminating in the 
central singularity at r = 0, the set /~(7) is a terminal indecomposable past (TIP) set 
as characterized by Geroch, Khronheimer and Penrose^^. The boundary of this past set 
is then generated by null geodesies which are endless in future, terminating at the central 
singularity, by the theorem (2.3) of Ref. [15]. It follows that infinitely many families of 
null geodesies terminate at the singularity in future. If the collapse develops from the initial 
spacelike surface t = 0, all these null geodesies meet this surface in the past in the TBL 
models. Consider now any event p in 1^(7). This later set being a TIP, the past of p, 
I~{p) must be contained in /~(7). Thus all the past directed nonspacelike curves from p 
must meet the initial surface t — Q, and so the set /~ (7) becomes a subset of the domain of 
dependence (see, e.g., Ref. [16] for definitions and further details) of the initial surface t — Q 
in the spacetime. Then, given any event p in and any other event q in the future of 

p in the same set, there exists a nonspacelike geodesic from p to q. Choosing now q to be 
closer and closer to the singularity, and in the limit approaching the same, we see that there 
are nonspacelike geodesies from p entering every neighbourhood of the central singularity. 
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and thus terminating in future only at the singularity. Thus, there will be infinitely many 
future directed nonspacelike geodesies, both null and timelike, which terminate at the central 
singularity in future. 

One can also analyze equations (10), (11) further to understand better the structure of 
timelike geodesies terminating at the central singularity in future. Another such solution 
describing ingoing timelike causal particles is given, in the neighbourhood of the central 
singularity (for n < 3), as 

R{t,r) = Xor^+2"/^ U' = 1 + ^ f :r-^V^'^'"^^ + ^'(r'"), (14) 



8 Vl -n/3 



where 



£1 
4 ' 

npi 

(3+n)(n-l)' 



n 



(2TJ2(3-n)), n^l, 



2_ ^Pi \ I Pc 



(3 + n)Xo='/VV3Xo' 

The trajectories described by the above equation actually represent a family in the spacetime 
of infinitely many ingoing radial timelike geodesies {9 — const, = const), terminating at 
the central singularity in future. In fact, the above timelike geodesies form a part of a larger 
family of timelike geodesic curves in the (t, r) plane that terminate in the singularity and 
are given in the near regions of the singularity by 



i?(t,r) =Xor^+'"/^ (16) 

where X^^^ — 2^^, C 7^ 1 is constant, and ± represent outgoing and ingoing geodesies. 
Figure 1 describes a typical radial timelike geodesic terminating at the central singularity. 
For the case n > 3 the singularity is already known to be strong^'^, and families of null and 
timelike geodesies terminate at the singularity. We have discussed the marginally bound 
(i.e., / = 0) case above for the sake of simplicity and clarity, however, for the nonmarginally 
bound case as well one gets similar results. As pointed out earlier, the solutions describing 
the particle trajectories in the {t, r) plane terminating in the singularity are not limited just 
to the solutions given by Eqs.(12)-(16), but there will be many more such families. 

The scenario which results from the above considerations is that a single ingoing timelike 
geodesic given by Eq. (12), and infinitely many other radial nonspacelike geodesies always 
terminate in future at the central singularity. This is regardless of the form of the initial 
data, which in the dust collapse is in the form of the initial distribution of the density and 
velocities at the onset of collapse. To understand better the structure of families (outgoing 
or ingoing) at the singularity, we note that in case p'(0) 7^ or p"(0) 7^ there are always 
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Figure 1: Typical ingoing radial timelike geodesies terminating at the central (r=0) sin- 
gularity in the TBL model. The analytic approximation is also shown. Here f — and 
p — 3 — 4.8r. (a) Different geodesies in the (r,R) plane, (b) C/*(r) for C—2 geodesic. 



families of ingoing and outgoing radial timelike geodesic curves that terminate at the naked 
central singularity. In cases where p'(0) = p"{0) = 0, p"'(0) ^ there is a critical value of 
p'" which determines if there would be outgoing null geodesic families; however, the ingoing 
family would still be there. Thus, in all cases (naked or covered), there is a nonzero measure 
of families of both timelike and null particles which terminate at the singularity. 

In order to access the strength of the singularity, consider the quantity ^ = RabU"'U^ 
along the central particle following the radial timelike geodesic 7(5) given by Eq. (12). The 
energy density throughout the spacetime is given by 

,U r) = ^ = /^(^)^°(^)^ (17) 

^'^ Rm' [t,{r)-t]-{[to{r)-t] + {2/3)[rt',{r)/to{r)]t} ^ ' 



Using this we get 



Since ^ diverges as the inverse square of the proper time, it follows that the timelike geodesic 
7(5) terminates in a strong curvature singularity. This shows, as per the criterion of strength 
of the singularity given above, that the central singularity forming in the dust collapse is 
always a strong curvature singularity. It follows, by continuity, that the curvature growth in 
the past of the singularity will be equally powerful in the limit of approach to the singularity 
along the other nonspacelike trajectories as well which meet the singularity with a tangent 
closer to that of the timelike geodesic 7(5) considered above. 
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A similar calculation as above for the particles following radial timelike trajectories given 
by Eq. (12) yields 



lim [{s - So)'*) = Jim ((s - sofRabW'U') 



2 



(19) 
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For the sake of completeness, we note that in the case of non-marginally bound dust collapse 
as well, where / 7^ 0, the quantity * along the timelike geodesies 7(s) given by Eqs. (12) 
and (14) behaves in the neighborhood of the singularity as 



As is wellknown, the stability criteria in general relativity have still not been satisfactorily 
formulated. However, certain norms for stability have come to be well accepted. Amongst 
them is the stability of the occurrence of phenomena against small perturbations in the 
symmetry of the spacetime. Recent work^^ on the perturbations of the TBL models shows 
the stability of the naked singularity in this sense. Another important stability criterion 
for any physical phenomena, which is the result of time development of initial physical 
configurations (e.g., the initial densities and velocities), is the stability against perturbations 
of these physical parameters. This is important from the point of view of observations and 
measurements, and also in the context of cosmic censorship, where the role of initial data 
leading to a naked singularity is quite important. 

To examine the stability of naked singularities from such a perspective, consider for 
example the marginally bound collapse scenario, in the case when for a given (at least 
differentiable) density profile p(r) the singularity is naked. Consider perturbations to this 
initial data of the form p ^ p + 6p. Since any such generic infinitesimal perturbation 
would require the introduction of nonvanishing first or second derivatives of the density 
[i.e. p'(0) 7^ or p"(0) 7^ 0; it is known^^ that for realistic stellar models p"(0) < 0], it 
follows that the singularity always continues to be naked. Furthermore, as shown above, 
in all these cases the strong curvature condition is always satisfied regardless of the initial 
data. It follows that the perturbations of initial data do not alter the occurrence of a strong 
curvature naked singularity, which is a stable phenomena in that the singularity remains 
not only naked but strong also. Similarly, in nonmarginally bound collapse also, a generic 
perturbation of the initial data consisting of the initial density and velocity profiles will leave 
the naked singularity to be strong. 

We thus see that the naked singularity developing in dust collapse is stable against such 
perturbations. What is seen is that the central singularity is always a strong curvature 
singularity, and infinitely many timelike geodesic trajectories terminate at this singularity, 
whether covered or naked. The important aspect here is the generality of this feature in- 
dependent of the initial density and velocity profiles. It is of course known from earlier 
studies that this singularity, especially when naked, does exhibit a directional behavior for 
the growth of curvature in the cases such as above, and also in other collapse models such 



RabU^U' oc 



1 



(20) 
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as the Vaidya spacetimes, along the radial null geodesies families terminating at the naked 
singularity, in the sense that the curvature growth may be different along different families. 
However, even in those cases the naked singularity is always strong as per the criterion given 
by Krolak^'^. Further, in a recent development^^, it has been shown that the redshift of rays 
from the naked singularity is finite or infinite depending upon how the density of the star 
decreases away from the center. Since the considerations here show the nature of the sin- 
gularity to be stable and strong, the central naked singularity formed in collapse becomes a 
good candidate for a physical phenomena likely to have an observational signature. Though 
we have only considered dust here, these results seem to hold for more general equations of 
state as well. This will be discussed elsewhere. 
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